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SEaiON  I 


INTRODUCTION 


The  investigation  of  hydrocarbon  exhaust  plume  characteristics  is  of 
interest  to  the  Air  Force  for  air-to-air  missiles  and  target  systems.  IR 
seeker  response  during  terminal  guidance  is  at  least  partially  determined 
by  the  spatial  distribution  of  plume  radiation. 

The  optical  and  IR  signature  of  the  plume  is  determined  by  the 
temperature  distribution  and  the  characteristic  frequencies  (spectrum 
lines)  of  the  radiating  species.  The  radiating  sjjecies  are  the  combustion 
products  and  radicals  formed  by  the  chemical  reactions  in  the  rocket 
plume.  For  this  reason,  experimental  determination  of  these  properties 
is  of  utmost  importance.  In  addition,  reliable  analytical  models  of 
hydrocarbon  rocket  or  jet  engine  plumes  need  to  be  developed  to  predict 
the  behavior  of  new  or  unknown  engine  designs. 

Optical  measurement  of  plume  properties  (temperature  and  intensi¬ 
ties)  is  mandatory  for  the  validation  of  these  models.  Another  way  to 
check  these  models  is  to  measure  the  concentration  of  the  combustion 
products  and  radicals  formed  in  the  rocket  plume  and  the  relative  changes 
in  concentration  of  these  species  if  the  oxygen  to  fuel  (0/F)  ratio  is 
changed.  This  information  can  also  be  obtained,  by  spectroscopic  measure¬ 
ments.  Consequently,  the  primary  objective  of  this  project  was  to  produce 
temperature  maps  of  the  exhaust  plume  generated  by  a  small  rocket  engine 
burning  RPl  with  gaseous  O2  as  the  oxidizer.  After  having  obtained  the 
temperature  distributions,  previous  measurements  of  spectral  intensities 
of  radicals  could  be  converted  into  relative  changes  in  concentration  as 
a  function  of  oxygen  to  fuel  ratio.  The  method  chosen  for  the  temperature 
measurements  was  the  line  reversal  technique  with  sodium  used  as  the 
emitting  species.  The  temperatures  obtained  with  such  a  technique  are 
most  accurate  when  obtained  from  a  homogeneous  plasma  of  constant  temper¬ 
ature,  density,  and  species  concentration.  When  such  conditions  are  not 
met  (such  as  in  rocket  exhaust  plumes)  the  measurements  produce  an  average 
temperature  due  to  the  existence  of  temperature  gradients.  Techniques 
are  proposed  for  unfolding  these  average  temperature  profiles  to  produce 
corrected  radial  temperature  distributions  of  the  plume.  The  investigations 
reported  here  are  the  first  experimental  application  of  these  techniques 
to  a  rocket  plume. 

Average  temperature  maps  were  made  of  plumes  produced  by  both  oxygen- 
rich  (0/F  =  4.5)  and  fuel -rich  (0/F  =  2)  mixture  ratios.  For  the  oxygen- 
rich  conditions  (0/F  =4.5)  both  radial  centerline  temperature  measurements 
were  made. 

A  corrected  (unfolded)  radial  temperature  distribution  was  generated 
for  one  plume  cross  section  at  the  0/F  =  4.5  condition.  The  unfolded 
temperature  profile  was  compared  with  the  average  temperature  profile 
obtained  at  the  same  location.  This  procedure  provides  an  estimate  of 
the  error  incurred  in  applying  standard  line  reversal  techniques  to 
non- homogeneous  rocket  plumes. 


1 


Due  to  unstable  burning  conditions  and  poor  plume  definition  at  the 
fuel-rich  mixture  ratio  (0/F  =2)  only  a  centerline  temperature  distribution 
was  produced.  However,  the  experience  gained  with  unfolding  of  the  0/F  = 
4.5  tranperature  profile  makes  it  possible  to  estimate  the  true  temperature 
of  these  measurements  even  without  undergoing  an  unfolding  procedure. 
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SECTION  II 


THE  MEASUREMENT  METHOD 


1.  LINE  REVERSAL  TEQINIQUE 

For  a  homogeneous  plasma  in  thermodynamic  equilibrium,  the  radiation 
field  is  given  by 


'x  - 


1 _ 

exp(hc^/AkT)-l 


(1) 


p 

where  (T)  is  the  Planck  function  and  T  is  the  plasma  temperature . 

In  case  of  ccanplete  thermodynamic  equilibrium,  the  plasma  is  a  blackbo- 
dy  ^d  its  radiation  is  only  a  function  of  surface  area  and  temperature. 

It  is  not  a  function  of  any  material  constants.  A  rocket  plume  is  not  such  a 
radiator,  it  is  only  local  in  thermodynmaic  equilibrium  (LTE).  In  this  case, 
the  radiation  is  determined  by  the  volume  of  the  plasma,  its  temperature  and 
a  materials  constant,  usually  expressed  as  an  absorption  coefficient.  Consi¬ 
der  a  homogeneous  slab  of  plasma  in  local  thermodynamic  equilibrium  (LTE)  as 
indicated  in  Fi^re  1.  At  any  given  point  inside  this  plasma  the  radiation 
field  is  determined  by  two  competing  factors.  The  plasma  emits  energy  of 
amount  ex  erg/ cm3  sec  per  unit  wavelength,  where  ex  is  the  emission  coeffi¬ 
cient.  It  also  absorbs  energy  at  the  rate  k(X)Ix(r)  where  k(X)  is  the  ab¬ 
sorption  coefficeint  of  the  plasma  at  a  particular  wavelength  and  Ix(r)  is 
the  radiation  intensity  of  some  position  r  in  the  plasma  slab..  The  rate  of 
change  of  intensity  through  the  slab  is  a  function  of  these  two  competing 
terms  and  is,  obviously,  given  by  their  difference 


dr 


-  k(X)l^(r) 


(2) 


It  is  assumed  that  the  plasma  slab  under  consideration  is  homogeneous,  which 
means  that  temperature,  emission  coefficient  and  abso^tion  coefficient  are 
not  a  function  of  r.  Of  course,  this  is  only  an  initial  assumption  for  dis¬ 
cussion  of  the  method.  In  a  real  plasma,  like  a  rocket  plume,  these  proper¬ 
ties  are  indeed  a  function  of  r.  In  this  case,  it  will  be  necessary  to  divide 
the  plume  into  rings  small  enough  that  it  can  be  assumed  that  within  one 
ring  the  plasma  is  homogeneous. 

The  assumption  of  LTE  requires  that  Kirchoff’s  Law  is  valid.  Kirchoff's 
Law  establishes  a  connection  between  and  k  (X).  It  states  that  the  ratio 
of  emission  coefficient  to  the  absorption  coefficient  is  only  a  function 
of  temperature  and  is  equal  to  the  Planck  function.  The  result  being  that 
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e.  =  k(A)Br(T).  C^xCT)  Planck  function 

of  the  plasma]  (3) 

Inserting  this  result  into  equation  (2)  yields 

dl,(r) 

-J-—  =  k(X)Bj(T)  -  k(X)l^(r) 

P 

=  k(X)tB^(T)  -  I;^(r)] 

This  may  be  integrated  utilizing  a  change  of  variable 

B^(T)  -  I^(r)  =  K 
3r 

(r)  =  B^  (T)  -  K 

md  *  ^ 

dlj^(r)  -  -dK 

Inserting  equation  (6)  and  equation  (8)  into  equation  (5)  we  have 

-  .  kWK 
^  =  -k(X)dr 

Integration  of  equation  (10)  produces 

—  =-k(X)/^dr 
o 


In  K  =  ■-k(X)r  +  In  K 


where 


K  =  Ce 


-k(X)r 


Utilizing  equation  (6)  yields 


I,  (r)  «  B,  (T)  -  Ce 


-k(X): 


Letting  ^  at  r  =  0,  C  may  be  evaluated 


=  ■  “x<«  -  'x,o 
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Finally,  inserting  equation  (16)  into  equation  (15)  yields 

-k(A)r 

l^{r)  =  B^(T)  -  [B^(T)  -  ®  (17) 

or 

P  -k(A)r 

I^{r)  =  B^(T)I  1  -  e  ®  (18) 

In  the  technique  of  line  reversal  temperature  measurement  a  standard 
illumination  source  is  imaged  upon  the  boundary  of  the  plasma  field.  The 
image  of  the  source  is  calibrated  against  a  standard  blackbody,  thereby 
giving  an  quivalent  brightness  temperature  for  the  illuminating  radiation 
yin  Figure  1.  The  incident  radiation  intensity  can  thereby  be  given  as 

h.o  - 

c 

where  B^(T)  is  the  Planck  function  for  the  source  image.  The  temperature 
of  the  standard  source  is  adjusted  until  the  output  intensity  equals  the 
incident  intensity  or 


I^(r)  =  B^(T)  I  1  -  *  B®(T)  (20) 


for  r  =  £  we  have 


B^(T)  = 


bJ(t) 


(  1  - 


e  J 


„S,„,  -k(X))l 
B^(T)e 


the  final  result  being  that 


B®(T)  =  B^(T) 


(21) 


(22) 


Since  the  Planck  function  is  a  function  only  of  temperature,  this  implies 
that  the  plasma  and  source  image  temperature  are  equal 

Tp  =  Tg  (23) 

The  above  technique  is  correct  only  for  homogeneous  slabs  in  local 
thermodynamic  equilibrium.  In  other  words,  the  states  of  atoms  or  molecules 
emitting  the  wavelength  region  which  is  used  for  the  line  reversal  techni¬ 
que  have  to  be  excited  according  to  a  Maxwell-Boltzmann  distribution. 

For  a  good  precision  of  the  measurement  it  is  also  required  that  the  absorp¬ 
tion  coefficient  is  of  a  certain  minimum  value.  Although  the  absorption 
coefficient  is  not  explicitly  measured,  the  method  depends  on  the  fact  that 
absorption  takes  place  in  the  plasma  [see  Equation  (2^  . 
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Absorption  requires  that  the  lower  state  of  the  spectrum  line  used  for 
the  method  is  appreciably  populated.  This  can  be  guaranteed  best  by  using 
as  the  lower  state  the  ground  state  of  the  atom.  For  this  reason  mostly 
resonance  lines  are  used  for  the  line  reversal  method.  Only  a  few  chemical 
elements  have  resonance  lines  in  the  visible  spectrum  region  (e.g.,  Na, 

KJ. 


2.  TEMPERATURE  UNFOLDING  PROCEDURE 

Temperatures  are  not  additive.  Therefore,  if  a  temperature  gradient 
exists,  regular  Abel  unfolding  techniques  which  are  based  on  the  fact  that 
intensities  are  additive  cannot  be  used. 

However,  the  quantity  t  =  k*8.,  Sometimes  referred  to  as  an  optical 
length,  is  additive.  A  temperature  unfolding  technique  can  be  based  on  this 
fact. 

Therefore,  a  radial  temperature  distribution  may  be  obtained  by 
dividing  the  plume  cross  section  into  thin  concentric  homogeneous  rings. 
Figure  2  illustrates  such  a  division  where  a  three-ring  section  is  utilized. 
Each  ring  is  assumed  to  contain  a  homogeneous  plasma  with  Planck  funtion 
Pn 

B,  (T)  and  absorption  coefficient  k  (X) .  This  assumption  requires  the 
A  n 

existence  of  LTE  with  Kirchoff's  Law  being  valid  in  a  narrow  spectral  range. 
The  accuracy  of  such  a  layering  technique  is,  of  course,  dependent  upon 
the  thickness  of  the  slabs,  therefore  A  is  kept  small, 

p 

The  calculation  of  the  b^''(t)  is  straightforward,  provided  the  outer¬ 
most  ring  temperature  is  known.  This  temperature  may  be  obtained  by  a  line 
reversal  temp)erature  measurement.  The  method  of  determining  the  remaining 

Pn 

(T)  is  now  explained. 

For  simplicity,  consider  a  plasma  divided  into  three  rings  as  seen 
in  Figure  3.  Consider  incident  radiation  of  intensity  L  .  The  total 

intensity  on  the  right  side  of  the  plasma  is  1.^  Let  Ip  "  be  the 

intensity  of  radiation  emitted  by  the  plasma  itself.  Then 


+  I,  e-Tr 
X,On  ^ 


(24) 


where  x  is  the  optical  thickness  of  the  plasma  for  a  given  path  and  is 
given  by^ 

n 


-  2  2:k  (X)Jln 
n  .  m  **, 


bifI 


»m 


(25) 


For  the  outermost  ray  (outside  ring  only  traversed)  we  have  that 

i_  =  i_  +  I. 

X,o 


(26) 


6 


vrtiere 


Tl  . 


or 


‘l'*’  ■  2i 


ln[ 


1.1 


-  1« 
T  P 
1  1 


1 


(27) 

(28) 


TTie  intensities  in  equation  (28)  are  easily  determined  and  ^  is  measure- 
able  from  the  plume  geometry.  Knowing  kj(X)  the  remaining  k^U)  "'ay  be 
determined. 


where 


2,2 


‘  K — C"' 


kj^(X)£2 

■  “2,2 


T2  -  2[k2(X)£2,2 


(29) 


(30) 


and 


wliere 


1C3(X)  - 


,  ^X.O,  £,  ,k  (X)  +  £  -k  (X) 

In  [r - = — -  I  - 0 -  ] 

3.3  ^T3  ^P3  S.3 


T3  .  2t£3^^k^(X)  +  il3^2’‘2<^)  ^  S,3’^3‘^>J 


(31) 


(32) 


P 

The  problem  now  remains  to  calculate  the  As  stated  earlier,  the 

temperature  and  thereby  the  Planck  function  of  the  outer  ring  are  deter¬ 
mined  through  line  reversal  measurements.  Therefore  Bfl(T)  is  known. 

A 

Consider  now  the  ray  that  passes  through  the  outer  two  rings  (Figure  4) . 
Inside  the  first  ring  the  intensity  may  be  written  as 

I(r^)  .  bJ1(T)  [  1  -  ^  “  ^2,1  (33) 

Inside  the  second  ring  we  have 

Kr^)  -  bJ2(t)  [  1  _  +  a  =  B  for  =  2l^  ^  (34) 


On  the  far  side  again  in  the  first  ring  we  have 
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-kx(A)r3 


(35) 


Kr^)  -  bJ1(T)  [  1  -  +  B 


At  reversal  for  ’t‘3  “  ^2  1 


.S2< 


'•'s’  -  'a,oj  ■ 


(36) 


where  Bx^(T)  is  the  Planck  function  for  the  source  image.  Since  the  objec¬ 
tive  of  this  technique  is  to  improve  upon  line  reversal  ten^jerature  measure¬ 
ments  the  B^(T)  are  known  fran  previously  generated  average  temperature 
profiles.  So  we  have  that 

B®2(t)  -  B^l(T)  (  1  - 
+  ^bJ2(t)  {1 

+  {bJ1(T)  [  1  -  e‘''l^^^^2,l] 

This  may  be  solved  for  B^2(t)  and  we  have 


bJ2 (t) 


Bf2(T)  -  B*’i(T)  II  - 
(  _ _  } 

e"’^2(^)*'2,l 


{bJ1(T)  [1  -  e-’^l<^’^2.1]  +  B®2(T)e’’'l<^>^2,l}  e-2’'2(^)«'2,2 

(1  -  e‘^’'2(^)^2,2] 


(38) 


Since  the  Planck  function  is  a  function  of  tenperature  only,  the  temperature 
of  the  second  ring  is  determined.  The  process  is  cascading  one,  i.e., 
each  consecutive  ring  temperature  requires  all  of  the  outer  ring  tenpera- 
tures.  For  the  path  through  three  rings  we  have  (Figure  5) 


Hr^)  -  B^i(T)  [  1  -  «  A  for  (39) 
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(40) 


Kr^)  «  B^2(t)  [  1  -  e'’'2(^)r2  j 

1^3)  =  B^3(T)  t  1  -  ] 

I(r^)  -  bJ2(t)  [  1  -  ] 

and  finally 

Kr^)  «  bJ1(T)  [  1  -  ] 

For  the  reversal  condition  we  have 


+  A  :=  B  for 

+  B  e  =  C  for  r3 

'+  C  e  «  D  for  r. 

4 

^  _  -k. (X)r_  „  - 

+  D  e  1  5  »  E  for  r^ 

-  B®3(T) 


^S,3  (41) 
S,2  C^2) 


(43) 


(44) 


and 


B^l(T)  [  1  -  +  D  -=  B®3(t)  (45) 

This  equation  contains  one  unknown,  namely  Bx^(T).  Thus  it  can  be  seen 
that  the  ring  temperatures  can  be  determined  by  the  above  procedure.  The 
technique  can  be  extended  to  any  number  of  rings  desired. 


3.  lETERMINATION  OF  CONCENTRATIONS 


As  pointed  out  in  the  introduction  the  measurement  of  concentrations 
of  radiating  species  provides  insight  into  the  chemistry  of  the  rocket  plume. 
In  case  of  atomic  species,  if  a  Maxwell-Boltzmann  distribution  of  states 
exists,  the  intensity  density  of  an  emission  line  of  a  spectrum  line  can 
be  related  to  the  particle  density  n^  (and  so  to  the  concentration)  by: 


I 


n  ^  hv 
o  U 


(46) 


vdiere 

g  »  statistical  weight  or  degeneracy  of  the  upper  state  of  the  line 
A  »  transition  probability 
hv*  energy  of  the  photon 
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U  ■  partition  function 
E  *  excitation  energy  of  the  line 
K  *  Boltzmann  factor 
T  *  ten^jerature. 

If  the  temperature  is  known,  measurement  of  the  intensity  will  produce  a 
value  for  n  ,  since  all  the  other  quantities  in  equation  (46)  are  known 
atomic  constants.  Measurement  of  relative  change  of  intensity  for  a  chang¬ 
ing  0/F  ratio  in  the  rocket  plume  can  be  interpreted  in  terms  of  the 
relative  change  in  concentration  of  a  certain  species,  again  under  the 
assumption  that  the  temperature  has  been  measured  by  some  means  (e.g.,  the 
line  reversal  method) .  The  complete  technique  is  described  in  Appendix  A. 

Most  of  the  radiating  species  in  the  rocket  plume  are,  however,  mole¬ 
cules  or  radicals.  In  this  case  a  modified  form  of  equation  (46)  which 
gives  the  intensity  of  a  line  in  a  rotational  band  must  be  used. 

I  «  (2J  +  l)e^.-I:Zltl2l 

0KT  (47) 


where 

J  *  rotational  Quantum  number 
2 

0  *  mr  (m  =  mass  of  molecule 

r  *  intermolecular  separation  distance) 
h  *  Planck’s  Constant 
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SECTION  III 
EXPERIMENTAL  SET-UP 


1.  ROCKET  ENGINE 

The  experiment  described  in  this  report  was  conducted  at  the  Range 
22  Stratochanfcer  Test  Facility  at  Eglin  Air  Force  Base.  The  capabilities 
of  this  installation  are  described  in  Reference  1.  The  experimental 
apparatus  consisted  of  a  small  kerosene/gaseous  oxygen  rocket  which  exhaust¬ 
ed  into  a  vacuum  chamber,  a  sodium  injection  system,  a  standard  illumina¬ 
tion  source,  an  imaging  optical  system,  and  a  spectroscopic  measurement 
system. 

The  rocket  engine  used  in  the  experiment  was  a  small  kerosene/gaseous 
oxygen  rocket  that  is  based  on  an  original  design  by  Astrosystems ,  Inc., 

Lake  Success,  New  York.  The  rocket  engine  injector  essentially  consists 
of  a  cylinder  with  a  0.055-inch-dian»ter  orifice.  Gaseous  oxygen  is 
introduced  into  the  injector  upstream  of  the  kerosene  and  at  a  higher 
pressure  than  the  kerosene.  This  design  permits  mixing  and  atomization 
to  occur  prior  to  injection  into  the  combustion  chamber.  In  the  wall  of 
the  combustion  chamber  is  a  port  for  injection  of  seed  material  (sodium 
chloride)  into  the  plume  for  application  of  the  line  reversal  technique. 
Typical  chamber  pressure  is  40  psi  requiring  that  the  seed  supply  line  be 
pressurized  at  a  higher  value  (80  to  100  psi) . 

The  engine  was  mounted  on  a  horizontal/vertical  translator  which  was 
electronically  controlled  outside  of  the  chamber  during  data  runs.  The 
translator  was  capable  of  engine  movements  of  1  mm  increments  along  both 
axes.  Associated  with  the  rocket  engine  was  the  control  panel  and  flow 
hardware.  Precision  micrometer  values  located  in  the  respective  lines 
permitted  the  fuel  and  oxygen  flow  rates  to  be  adjusted  to  achieve  a 
desired  oxygen  to  fuel  ratio  during  operation  of  the  engine.  The  combus¬ 
tible  mixture  is  initially  ignited  by  a  spark  plug.  Total  mass  flow  for 
the  engine  was  maintained  at  5  grams  per  second.  For  the  0/F  *  2 
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condition,  the  flow  rates  were  oxygen  3.33  grams  per  second,  RPl  1.67  grams 
per  second  and  for  0/F  =  4.5,  oxygen  4.10  grams  per  second,  RPI  90  grams 
per  second. 

A  sodium  chloride/water  solution  of  47  grams  per  liter  was  used  as 
the  plume  seed  material.  The  solution  was  force-fed  to  the  engine  through 
a  0.25  inch  copper  line  from  a  nitrogen  pressurized  steel  tank.  Tank 
pressure  ranged  from  80  to  100  psi.  Metering  was  accomplished  by  a  pre¬ 
cision  micrometer  needle  valve.  To  facilitate  metering,  a  0.018-inch- 
diameter  orifice  was  placed  ahead  of  the  combustion  chamber  port.  The 
orifice  also  acted  as  an  atomizer  for  the  sodium  chloride  solution  and 
prevented  extinction  of  the  plume  by  excess  water.  Consistent  sodium 
chloride  injection  was  insured  by  constant  monitoring  of  combustion  cham¬ 
ber  pressure  and  maintenance  of  this  pressure  at  approximately  40  psi. 

2.  OPTICAL  ARRANGEMENT 

The  objectives  of  the  optical  imaging  system  were  to  (1)  image  a 
tungsten  ribbon  lamp  filament  upon  the  rocket  plume  centerline  with  a 
magnification  ratio  equal  to  unity,  and  (2)  image  the  plume  and  filament 
upon  the  entrance  slit  of  the  spectrometer  with  a  convenient  magnification 
ratio.  A  schematic  diagram  of  the  imaging  system  is  shown  in  Figure  6. 

The  tungsten  filament  was  imaged  by  a  telephoto  astigmat  lens  with 
an  effective  focal  length  of  20  inches  and  f-number  of  5.6.  Unity  magni¬ 
fication  required  that  the  filament  be  located  a  distance  of  (2)  x 
(effective  focal  length)  from  the  principal  plane  of  the  lens.  Verification 
of  correct  distances  is  accomplished  by  accurately  measuring  the  filament 
size  and  adjusting  the  optics  to  obtain  unity  magnification. 

The  filament  image  and  plume  were  imaged  onto  the  spectrometer  entrance 
slit  by  a  12-inch-diameter  parabolic  mirror  of  focal  length  2.51  meters. 

This  is  facilitated  by  folding  the  beam  with  a  12 -inch-diameter  plane  front 
surface  mirror.  The  plume  and  filament  image  are  slightly  magnified  by 
this  system  with  a  magnification  ratio  of  1.05.  The  spectrometer  used  was 
a  0.75  meter  Czemy-Turner  mount  (manufacturer:  Jarrell-Ash) . 
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The  spectrometer  was  equipped  with  a  148- groove /mm  grating  blazed  for  5 
microns . 

A  500-element  silicon  multichannel  detector  unit  was  mounted  at  the 
exit  plane  of  the  spectrometer  and  the  grating  rotated  to  allow  illumination 
of  the  detector  array  by  the  sodium  doublet  (approximately  5890X) .  The 
multichannel  detector  unit  was  interfaced  with  an  Optical  Nfeiltichannel 
Analyzer  (CMA).  The  OMA  (manufacturer:  Princeton  Applied  Research)  allowed 
real  time  display  of  relative  incident  intensities  upon  the  silicon 
detector  array  via  oscilloscope  interfacing,  as  indicated  in  Figure  7. 

Such  relative  intensity  displays  could  be  stored  in  the  OMA  memory 
unit  and  plotted  for  permanent  record  with  an  X-Y  plotter. 

Figure  8  shows  a  typical  output  of  the  X-Y  plotter  for  pre-reversal 
and  post- reversal  conditions.  Figure  9  shows  how  actual  reversal  was 
established.  Once  a  reversal  condition  was  established  (center  part  of 
Figure  9)  the  lairp  ten5)erature  was  increased  slightly  (upper  part  of  Figure 
9)  and  decreased  slightly  (lower  part  of  Figure  9)  until  a  distinct  signal 
could  be  detected.  In  the  case  of  Figure  9,  this  required  an  increase  or 
decrease  of  the  shunt  voltage  drop  by  0.04  mV.  This  translates  into  change 
in  temperature  of  The  sensitivity  of  the  method,  which  is  the  rela¬ 

tive  error,  is  therefore  better  that  this  value. 

The  standard  light  source  consisted  of  a  calibrated  tungsten  ribbon 
lamp  with  its  associated  variable  high  current  dc  power  supply.  A  preci¬ 
sion  shunt  was  inserted  in  series  with  the  lamp  filament  and  the  voltage 
drop  across  the  shunt  was  monitored  by  a  dc  micro -voltmeter.  This  allowed 
calibration  of  the  lamp  with  an  optical  pyrometer,  the  result  being 
establishment  of  a  relationship  between  shunt  voltage  and  filament  bright¬ 
ness  temperature.  The  calibration  curve  is  shown  in  Figure  10. 

As  pointed  out  above,  in  the  line  reversal  technique  the  tungsten  rib¬ 
bon  lamp  filament  is  imaged  onto  the  boundary  of  the  plasma  field  produced 
by  the  rocket  plume.  The  filament  image  thereby  becomes  the  new  standard 
source  for  experimental  measurements.  This  necessitates  calibration  of 
the  filament  image  with  an  optical  pyrometer. 


13 


The  resultant  calibration  curve  vrtiich  relates  the  shunt  voltage  of  the 
filament  of  the  standard  lamp  to  the  brightness  temperature  of  the  fila¬ 
ment  image  at  the  centerline  of  the  plume  is  shown  in  Figure  11.  This 
calibration  curve  is  slightly  different  from  Figure  10,  due  to  the  fact 
that  the  standard  lamp  radiation  experiences  a  slight  absorption  at  window 
No.  1  when  entering  the  vacuum  chamber.  Of  course,  there  is  also  some 
absorption  in  window  No.  2.  But,  since  both  the  light  coming  from  the 
standard  lamp  and  light  coming  from  the  plume  are  absorbed  in  the  same 
way,  the  influence  of  window  No.  2  cancels  and  the  temperature  measurement 
is  not  affected. 
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SECTION  IV 


EXPERIMENTAL  RESULTS 

1.  TEMPERATURE  MEASUREMENTS 

A  photograph  o£  the  exhaust  plume  is  shown  in  Figure  12.  There  is  a 
visible  shock  structure  which  corresponds  quite  will  to  the  structure  ob¬ 
served  in  the  temperature  scans. 

These  results  are  shown  in  Figures  13  and  14.  The  spatial  resolution 

of  the  measurement  setup  was  1  mm.  Figure  13  is  a  plot  of  the  centerline 

temperature  of  the  plume  for  an  0/F  =  4.5,  which  is  a  lean  mixture.  The 

nozzle  exit  temperature  was  found  to  be  T^^  *  2503°K,  which  compares  very 

0 

well  with  the  theoretical  prediction  of  T  *  2570  K  reported  by  Vamos 

(Reference  2) .  The  temperature  rises  in  the  front  of  each  shock  drastically 

(from  1825°K  to  2245°K  within  1  mm  at  the  first  shock  front)  and  decreases 

0 

to  a  temperature  of  T  *  1827  K  which  is  fairly  uniform  between  all  shocks. 
Figure  14  shows  the  same  measurement  for  an  0/F  *  2.0  (rich  mixture). 

The  temperatures  are,  in  general,  lower  and  the  temperature  rise  in  the 
shock  front  is  not  nearly  as  pronounced  as  in  the  0/F  *  4.5  case.  These 
findings  have  to  be  compared  to  Figures  18,  19,  and  20,  which  indicate 
increased  formation  of  radicals  with  increasing  richness  of  expense  the 

fuel  mixture.  This  means  formation  and  excitation  of  the  radicals  at  the 
expense  of  randan  energy  of  the  plume  proper. 

The  line  reversal  measurements  indicate  that  the  plume  exhibits  a 
moderate  ten?)erature  gradient  in  the  radial  direction.  A  reversal  tem¬ 
perature  profile  for  one  cross  section  at  the  0/F  *  4.5  condition  is 
shown  in  Figure  15.  This  profile  was  unfolded  to  produce  a  true  radial 
tenperature  distribution  which  is  shown  in  Figure  16.  A  full  discussion 
of  the  unfolding  technique  including  actual  data  and  a  computer  program 
for  numerical  evaluation  of  the  data  are  contained  in  Appendix  B.  A 
large  number  of  cross  sections  were  scanned  and  this  material  is  compiled 
in  Appendix  C.  From  this  material  a  temperature  map  of  the  plume  was 
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constructed  for  the  case  of  0/F  »  4.5.  This  map  is  shown  in  Figure  17. 
Since  the  radial  temperature  gradients  were  moderate,  the  correction  ob¬ 
tained  by  the  unfolding  process  was  less  than  30  percent.  However,  the 
precision  of  the  line  reversal  technique  for  the  temperature  range  of  in¬ 
terest  for  this  investigation  is  +1  percent.  Therefore,  the  corrections 
resulting  from  this  unfolding  process  seem  very  worthwhile. 

2.  RELATIVE  CONCENTRATIONS 

The  visible  and  ultraviolet  spectra  were  obtained  in  the  rocket  plume 
previously  (Reference  3) .  The  observed  spectra  were  those  of  the  radicals 
OH,  CH,  and  C2  (Swan  and  higher  pressure  bandhead  at  5165^).  The  relative 
intensity  of  these  bands  was  measured  as  a  function  of  oxygen  to  fuel 
mixture  ratio.  In  all  cases,  an  extremely  sharp  increase  in  band  intensi¬ 
ty  occurs  vdien  the  mixture  ratio  is  reduced  below  0/F  =  3.5  (stoichiometric) 
and  reaching  a  maximum  at  0/F  »  1.8.  The  same  behavior  is  found  for  the 
total  CO2  in  the  4- to  5 -micrometer  infrared  emission  band.  These  data  show 
that  the  plume  radiation  is  primarily  due  to  external  combustion  between 
fuel -rich  exhaust  products  and  air. 

The  obtained  relative  intensities  were  converted  into  relative  concen¬ 
trations  using  the  reported  temperatures.  Details  of  the  conversion 
technique  are  contained  in  Appendix  A.  The  results  based  on  the  above 
three  bandhead  wavelengths  are  shown  in  Figures  18,  19  and  20.  In  all 
cases,  an  extremely  sharp  increase  in  band  radiation  intensity  occurs  as 
the  0/F  ratio  is  reduced  below  3.5  (stoichiometry).  This  variation  with 
0/F  ratio  is  similar  to  that  observed  in  molecular  CO2  emission  in  the 
infrared  band  between  4  and  5  microns,  although  the  rate  of  increase  in 
visible  band  radiation  is  much  larger.  The  previously  measured  (Reference 
3)  infrared  radiant  intensity  variation  with  mixture  ratio  is  reproduced 
in  Figure  21  to  show  the  con^jariswi. 

The  presence  of  strong  C2  and  CH  emission  are  indicators  of  fuel-rich 
flame  combustion.  Their  concentrations  begin  to  fall  off  below  a  mixture 
ratio  of  1.8,  as  does  the  infrared  radiant  intensity,  indication  that  the 
threshold  for  carbon  luminosity  may  have  been  reached.  However,  although 
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carbon  soot  is  visually  evident  in  the  exhaust  of  engines  under  these  fuel- 
rich  conditions,  it  is  apparent  that  the  amount  is  nonetheless  not  high 
enough  to  contribute  significant  infrared  radiation. 

The  data  thus  shov/  the  rocket  engine  radiation  to  be  primarily  due  to 
external  ccxnbustion.  An  increase  in  infrared  radiation  accanpanies  the 
increase  in  visible  radiation,  i.e. ,  the  increase  in  external  combustion. 

No  graybody  radiation  is  present  in  the  infrared  spectrum  under  these 
conditions  (Reference  3)  and  only  a  limited  visible  continuum  was  ob¬ 
served  over  the  shock  regions. 

OH  radical  was  determined  to  be  the  major  plime  emitter  at  visible  and 
ultraviolet  wavelenghts  at  high  mixture  ratios  and  beyond  the  region  of 
the  shocks.  None  of  the  major  molecular  coni)ustion  products,  such  as  H2O, 
some  of  which  are  primary  infrared  radii  tors,  were  observed 
to  emit  in  the  visible  or  ultraviolet.  The  0/F  ratio,  as  in  the  infrared 
wavelengths,  is  the  predominant  parameter  affecting  the  intensity  and  spec¬ 
trum  of  OH  radiation  in  the  UV. 
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Figure  3.  Three  Ring  Division  System 
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Figure  4,  Ray  Path  Through  Two  Rings 


Experimental  Setup  for  Sodium 

Line  Reversal  Plume  Temperature  Mapping 
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Figure  6.  Imaging  System 
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Figure  7.  Oscilloscope  Display  for  Sodiim  Doublet 
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Figure  9.  Sample  Data  Set 
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Figure  13.  Centerline  Tespersture  Distribution  (0/F  *  4.5) 
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Figure  14.  Centerline  Ten^rature  Distribution  (0/F  =2.0) 
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Figure  16.  Unfolded  Plume  Temperature  Distribution 


Figure  17.  Tenperature  Map  of  Plixae 


2  3  4  5  6 

Figure  18.  Relative  OH  Concentration  Versus  0/F  Ratio 
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Figure  19.  Relative  CH  Concentration  Versus  0/F  Ratio 
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Figure  20.  Relative  C2  Concentration  Versus  0/F  Ratio 
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OONCEmATIGN 


For  the  band-head  line  of  a  molecular  band  we  have  that 
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TTius,  the  connection  between  the  relative  intensity  and  relative  concen¬ 
tration  is  established. 

The  relative  intensity  distribution  is  obtainable  from  previous  mea¬ 
surements  (Reference  3) .  The  aie  tabulated  in  the  Appendix,  of  Reference 
4. 

A  distribution  of  the  plume  temperatures  for  0/F  -  2  through  0/F  ■  6 
is  obtained  by  fitting  a  straight  line  through  two  known  temperatures. 

These  taiqjeratures  were  taken  to  be  those  measured  by  sodium  line  rever¬ 
sal  at  the  first  shock  front  of  the  plume  at  the  0/F  -  2  and  0/F  ■4.5 
conditions. 
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APPENDIX  B 


PLUME  TEMPERATURE 
UNFOLDING  TECHNIQUE 
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UNFOLDING  TECHNIQUE 


The  plume  geometry  is  shown  in  Figure  B-1.  The  plune  cutoff  limit 
was  determined  photographically  to  be  7.26  mm  from  the  plume  centerline 
axis.  With  the  exception  of  the  outer  ring,  all  rings  have  a  thickness 
of  1  mm  and  data  paths  are  taken  through  the  center  of  each  ring.  Poor 
signal  to  noise  ratio  near  the  outer  boundary  of  the  plume  prevented  the 
use  of  a  seven-ring  system. 

All  path  lengths  were  determined  graphically  from  a  scaled  drawing 
of  the  plume. 

For  each  path,  the  following  procedure  was  followed:  (1)  intensity  of 
the  plasma  only  (Ip)  was  measured  (5889^)  and  (2)  the  tungsten  ribbon 
lan^  was  increased  past  the  reversal  point  and  the  baseline  (I^)  and 
minimun  (I^)  intensities  were  measured.  The  intensity  displays  and  the 
tabulated  intensities  appear  in  Figure  B-2. 

A  con^juter  program  was  written  in  order  to  simplify  the  analysis  of 
the  data  and  this  appears  at  the  end  of  this  appendix. 

The  outer  ring  is  taken  as  ring  1.  For  data  input 

Ip(n)  >S(N) 

I^(n)  -*-T(n) 

I^Cn)  -^E(n) 

The  TLAMP  (n)  value  corresponds  to  the  measured  line  reversal  tempera¬ 
tures  through  each  respective  path. 

The  computer  output  includes  a  list  of  the  ring  absorption  coefficients 
[kappa  (n)]  and  the  ring  temperatures  [TEMP  (n)]. 

For  the  plume  data  tabulated  here  the  KAPPA  (n)  are  seen  to  go  negative 
in  the  fifth  and  sixth  ring.  This  is  caused  by  the  fact  that  the  plasma 
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is  optically  very  thin  (K»t<l)  and  the  fact  that  errors  are  compounded  as 
the  unfolding  process  advances  toward  the  center  rmg.  Ring  temperatures 
may  be  obtained  only  for  positive  KAPPAS  so  the  generated  ring  temperatures 
are  useful  only  through  the  fourth  ring.  This  indicates  that  the  center 
temperature  of  the  plume  is  at  least  that  attained  in  the  fourth  ring 
(2697°K) . 
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/list 

0000  1  DIWCION  PATHC6,6),KAPPAC6)^EC6)^TC6),SC6),TERMC6)^BPLASC6),BLM«> 

0001  1C6),AC6),VC6),T0*»C6),TUW6) 

0002  2  PSlf^O 


ml 


MK>^iAior>eoo^pp^CM 


CM 


CM 


K\^iA«Arvoo9to*-icMK\^iA«or^eoo«Of-icMK\tftn<pr>«ooAeiMCMro 

OOOeOOO<-«f^f-if-4<Hi-4f-4iMtHiHeMCMMeMCMCMMeMCMMK\K>KkK\ 

ooooooooooooooooopooooooooooooo 

ooooooooooooooooooooooooooooooo 
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z «/»  — 

III 

y>  M 

sSs 


>  ^  lA  U>  I 
I  fA  KM 


> 

<  - 


I  tH  W  tA 
J"  ^  Jf 


5  S 

«  5  W  •  ^ 

^  ^z^'-i^^zt?' 

fknvpr<%apo^oWc<4fA. 
.^54f?d’>AlAIAtAI 


Pi 

III 


J-  IA(A 
lA  lA  tA 


tv.  00  CT» 

lA  lA  lA  ( 


^  § 


Pzfi^i 

/A  < 

.  lU  1-4  UJ  to  I 

!j||!j||; 

rA  .rf  tA  vO  I 
to  tO  to  to  < 


^IAtOr»OOOtQ<-)MtA4’IAtOr>«OOq>Of-)C><»A^IAtOr^OOOtOi-ICMrAd'IAlOA«OOOtO 

rArAK\IAtA(A^O'>^J.^J^^3'tAIAtAIAIAIAIAIAIAIAlOtOtOtOtOtOtOtOtOtOr» 

oooooooeooooeooooooooooooo  o  o  ooooeoooo 


K\rAK\IAtA»A4’J-^J‘J.* 

oooooooeoooo 

oooooooooooo 


ooooooooooooooooooo 


o  o  o  o 
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0071  69  OQNTINUE 

0072  70  (>073J=1,N 

0073  71  WAVE=CU»MB)«(l.E-8) 

0074  72  BUMPCJ>=Cl.l91E-5y«CWAVE»^5)/CEXPCl.4388/WAVE/TU^CJ))) 

0075  73  CONTINUE 


^  *H  T  /*> 
—  ?  1  Vi-3 

Z  <  "5^1 

In  1 

00  N.  -3  V./  I 


wrs  V./  XV 

M  ^  F"4 

I  ^  I  ^  £ 

V./  «-•  vf*  “3  m  w 

II  u  h 


Z  X  I  _L  ^ 

«  "3  V 

J  ^  UJ 

o5  I  «  >  /->  5 

XgSXii 


II  v 

WN  Z 


IP 

ssPs 


tnvo 

hN  r>. 


>  o  •-•  N 

.00  00  00 


0*0  >-t 
CO  0*0* 


\or^ooo*Or^o^K\0■*n\Dt^ooo*o<-^c^t<\■0>n^or^ooq* 

r^t^t^r>,aooooocococooooocono*o*o*o*^o*o*o*o*0* 

ooQooooooooeeeoooeoooooo 

oooooooooooooooooooooooo 
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0100  92  MESSAGE* EQUALS ' 

0101  93  WRITEC6)T0«»Cd) 

0102  94  CONTINUE 

0103  95  BO 

END  OF  WORK  FILE 


DATA  INPUT  CPATH  VALUES) 
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00 
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AMBOA  DESIRED  IN  >^STROM5 
60  REM>C5)LAMB 
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68  READC5)Tl/MPCJ) 


OUTPUT  CKAPPA  VALUES) 
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I 

i 

! 
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DATA  INPUT  CINTENSITY  VALUES) 
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034  32  READ(5)TCI)  040  38  READ(5)SCI) 


(Xn-PUT  CRING  TEMPERATURES) 
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Ray  1  X _ _ I  ^11=  4.83nm 
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Figure  B-1.  Plume  Geometry 


All  Intensities  in  Relative  IMits 


yi)  .  17 

I  (2)  »  17 
P 

Ip(3)  -  23 
Ip(4)  »  52 
I  (5)  =  66 
Ip(6)  =  70 


Plasma  ftily 
(No  Lan^) 


Il(1) 

-  195 

IjCl) 

m 

139 

lLt2) 

=  201 

Iy(2) 

s 

127 

1^(3) 

=  213 

1^(3) 

s 

153 

Il(4) 

=  396 

1^.(4) 

s 

310 

Il(5) 

=  679 

1^.(5) 

= 

660 

Il(6) 

=  881 

1^.(6) 

= 

860 

I 


I 


Lamp  Driven  Past 
Reversal  Point 


Figure  B-2.  Data  Collection  Technique 
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APPENDIX  C 

RADIAL  TEMPERATURE 
SCANS 
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(  Nq)  ajruBjaduisi  xesadAd^ 
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Figure  C-1.  Radial  Ten|)erature  Scan  (1.0  on  Aft  of  Rocket  Nozzle) 
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Figure  C-2.  Radial  Tenfjerature  Scan  (1.9  cm  Aft  of  Rocket  Nozzle) 
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Figure  C-3.  Radial  Temperature  Scan  (2.0  cm  Aft  of  Rocket  Nozzle) 
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Figure  C-5.  Radial  Teiqierature  Scan  (4  cm  Aft  of  Rocket  Nozzle) 
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Figure  C-6.  Radial  Tenqjerature  Scan  (5.0  an  Aft  of  Rocket  Nozzle) 


t)[^l  oanjeaaduiax  lesaoAoy 


64 


Figure  C-7.  Radial  Ten?)erature  Scan  (7.5  an  Aft  of  Rocket  Nozzle) 
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Figure  C-8.  Radial  Ten^rature  Scan  (8.5  cm  Aft  of  Rocket  Nozzle) 
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Olp)  aan:iBaadui9i  iBsasAay 
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Figure  C-11.  Radial  Tenperature  Scan  (15  cm  Aft  of  Rocket  Nozzle) 
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Figure  C-12.  Radial  Temperature  Scan  (20  cm  Aft  of  Rocket  Nozzle) 
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Figure  C-13.  Radial  Tenqjerature  Scan  (25  an  Aft  of  Rocket  Nozzle) 


Figure  C-14.  Radial  Temperature  Scan  (30.0  cm  Aft  of  Rocket  Nozzle) 


INITIAL  DISTRIBUTION 


HQ,  USAF/SAMI 
HQ  USAF/RDPA 
HQ  USAF/RDQRM 
HQ  USAF/XOOFA 
AFSC/INA 
AFSC/SDA 
AFSC/DLCAW 
AFIS/IOT 
AFSC/DPSL 
AFSC/SDWM 
TAC/DRA 
TAC/XPSY 
AFAL/AA 
AFAL/TEM 
AFAL/RWM 
AFAL/RW 
ASD/YHEV 
ASD/XRG 

ASD/YFEI  (ARMAMENT) 

ASD/ENO 
ASD/ENYW 
ASD/ENASA 
ASD/YPEX 
ASD/RWS 
ASD/RWR 
ASD/ENFEA 
ASD/SD  (TECH  DIR) 

ASD/SD7 
ASD/SD5EI 
ASD/SD4T 
AFFDL/FE 
AFFDL/FGL 
AFFDL/FX 
AFFDL/FY 
AFML/MX 
AFML/LL 
AFML/MBy 
AFML/I 
AFLp/mVM 
OO-ALC/NTWMP 
6585  TG/GDP 
AUL/LSE  71-249 
ATC/XPQS 

NAV  AIR  SYSTEMS  COMD/AIR-5323 
NAV  AIR  SYSTEMS  COMD/AIR-5324 


ODDR8E/TST5E 
DARPA/TIO 
NAVAIR  SYS  COMD 
USAF/AFSC  LIASON  OFFICE 
C^WANDER  NWC/CODE  456 

CX»MANDER  NWC/GODE  533 

CORIANDER  NWC/GODE  4063  y 

(BWANDER  NWC/CODE  40903X 

COLDER  NWC/CODE  330. 

CONMANDER  NWC/CODE  3, 

REDSTONE  SCI  INFO 
DSN  WEAPONS  LAB 
CG  USAMICOM 
DDC-TC 

CINCPACAFZ/GFW 
ADTC/PP 
ADTC/' 

ADTi 
?/DR 
tfC/ERW 
'TAWC/TX 
AFATL/DL 
AFATL/DLB 
,AFATL/DLMI 
FATL/DLNM 
FATL/DLOU 
FATL/DLY 
IfATL/DLYA 
ATL/DLYW 

afAtl/du 

Aim/DUF 
AFA’I\/DUC 
AFATIYDUA 
AFATL/^LODL 
ADTC/S 
ADTC/S 
ADTC/SDE 
Ain’C/SD7 
ASD/SDO 
ASD/SDM 
AFFDL/FGC 
AFML/LPO 
USAF  ACADEMY 
SAC/XHHN 
SAC/DOXT 
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SAC/LGWC  1 
AFAL/NVA-679A  2 
ASD/SD-65  1 
TAWC/TXA  1 
PHILO)- FORD  1 
ROCKWELL  INTERNATIONAL  1 
HUGHES  AIRCRAFT  00.  1 
MARTIN-MARIETTA  1 
RAND  CORPORATION  1 
AFRPL/MK  1 
GENERAL  DYNAMICS  1 
TAC/INAT  1 
ASD/XRP  1 
HQ  USAFE/DOQ  1 
USA  TRADOC  SYS  ANALYS  ACTVY  1 
HQ  PACAF/DOOFQ  3 
OOMIPAC/I-232  1 
ASD/ENESS  1 
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